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SUMMARY 

thp t^hnioue of differential display reverse iranscription-polymerase chain reaction (ddRT-PCR) has been used to produce unique 
«mffl« of Vreeulated and down-regulated gene expression in the liver of male Wistar rats following short term exposure to the 
™ r/notoxic hepatocarcinogens, phenobarbital and WY-14,643. Animals were treated for 3 days, whereupon their livers were 
SSTS snlpTozen. niRNA was prepared from the livers and used for ddRT-PCR. Individual bands from the different,* 
ojs^y7were extracted and cloned False positives were eliminated by dotblot screening and true positives then sequenced and 

identified. 



INTRODUCTION 

Safety evaluation of. new chemicals usually necessi- 
tates the examination of genotoxic and carcinogenic 
potential using short-term in vitro and in vivo geno- 
toxicity assays augmented by chronic bioassay tests. 
The short-term assays have proved useful in the early 
identification of potential genotoxic carcinogens, but 
their value is limited by observations which suggest 
that approximately 60% of chemicals identified as car- 
cinogens in life-exposure studies produce mainly 
negative findings in short-term genotoxicity tests (1,2). 
Thus, there is currently no reliable and rapid means of 
evaluating the carcinogenic risk of new chemicals 
which fall into this latter group of compounds, termed 
non-genotoxic (or epigenetic) carcinogens. 



Please send reprint requests to : Dr John Rockett. Molecular 
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It is now evident that non-genotoxic carcinogens 
constitute a group of chemicals which axe not only di- 
vergent in their interspecies toxicity, but also demon- 
strate different target organ selectivities and mecha- 
nisms of action (3,4). Elucidation of the molecular 
mechanisms underlying non-genotoxic carcinogenesis 
is currently underway, but the picture is still far from 
complete. It is anticipated that a better understanding 
of the early changes in genetic expression following 
exposure to non-genotoxic carcinogens will aid devel- 
opment of experimental strategies to identify cellular 
markers which are diagnostic for this type of toxicity. 

Subtractive ddRT-PCR is a recently developed 
technique which facilitates the preferential amplifica- 
tion of gene products that demonstrate altered expres- 
sion in target tissue(s) following exposure to chemical 
stimuli. Furthermore, using this technique, no prior 
knowledge of the specific genes which are up/down 
regulated is required. In the current study, we have un- 
dertaken to develop a specific and rapid assay for non- 
genotoxic carcinogens using the technique of ddRT- 
PCR. This has allowed us to identify characteristic 
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patterns of gene regulation following administration of 
two different non-genotoxic carcinogens (phenobarbi- 
tal and Wy- 14,643) and the subsequent identification 
of individual gene species which are regulated by this 
xenobiotic treatment. 

MATERIALS AND METHODS 
Animals and treatment 

Phenobarbital (BDH, Poole, UK; 100 mg/kg/day) or 
[4<hloro-6-(2 1 3-xylidino)-2-pyrimidinylthio] acetic 
acid (Wy-14,643) (Campo, Emmerich; 250 
mg/kg/day) was administered by gavage to groups of 
3 male Wistar rats (150-200 g) on three consecutive 
days, whilst control animals received nothing. All ani- 
mals had free access to food (rat and mouse standard 
diet, B&K Universal, Hull, UK) and water. The ani- 
mals were killed on the fourth day, whereupon their 
livers were excised, sliced into 0.5 cm cubes, snap fro- 
zen in liquid nitrogen and then stored at -70°C. 

mRNA extraction 

Up to 0.25 g of each frozen liver sample was ground 
under liquid nitrogen using a mortar and pestle. 
mRNA was extracted from the ground liver using 
Promega's PolyATtract® System 1000 (Promega, 
Madison, WI, USA) according to the technical man- 
ual. The mRNA was DNase-treated (Promega, final 
concentration 10 U/ml) before phenol/chloroform ex- 
traction and ethanol precipitation. The mRNA was re- 
suspended at a final concentration 500-1000 ng/uJ. 

ddRT-PCR 

This was carried out using the PCR-Select™ cDNA 
Subtraction Kit (Clontech, Palo Alto, CA, USA) ac- 
cording to the manufacturer's instructions. Final PCR 
reactions were run on a 2% Metaphor agarose (FMC, 
Rockland, MD, USA) gel containing ethidium bro- 
mide (Sigma, Dorset, UK) and then overstated for 30 
min with SYBR Green I DNA stain (FMC, 1:10 000 
dilution in TAE). 

Band extraction and cloning 

Each discernible band from the differential display 
pattern was extracted from the gel with a scalpel and 



the DNA eluted using a Genelute™ Agarose Spin Col- 
umn (Supelco, Bellefonte). An aliquot of the eluted 
DNA (5 u.1) was re-amplified using the original ddRT- 
PCR nested primers and electrophoresed on a 2% 
agarose gel. The re-amplified band was extracted from 
the gel (as above) and the eluted DNA ligated directly 
into the TOPO TA Cloning® vector (Invitrogen, 
Carlsbad) before transformation in Escherichia coli 
TOP10F One Shot™ cells (Invitrogen). 

Stage 1 screening 

Twelve transformed (white) colonies from each band 
were grown up for 6 h in 200 LB broth containing 
ampicillin (Sigma, 50 u.g/ml) and 1 |il of this ampli- 
fied by PCR reaction (as specified in ddRT-PCR tech- 
nical manual). One quarter of the completed reaction 
was electrophoresed on a standard 2% agarose gel and 
one quarter on a 2% agarose gel containing HA Yel- 
low (Hanse Analytik GmbH, Bremen, Germany, 1 
U/|il) to discern the different cloning products. The re- 
mainder was used to prepare duplicate dotblots on Hy- 
bond N+ (nylon) membranes (Amersham, Little Chal- 
font, UK). Cultures containing different cloning prod- 
ucts were grown up and a plasmid miniprep prepared 
from each (Wizard Plus SV Minipreps DNA Purifica- 
tion System, Promega) according to the manufac- 
turer's instructions. 



Stage II screening 

The duplicate dotblots were probed with: (a) the final 
differential display reaction; and (b) the 'reverse-sub- 
tracted* differential display reaction. To make the 're- 
verse-subtracted' probe, the subtractive hybridisation 
step of the ddRT-PCR procedure was carried out using 
the original tester cDNA as a driver and the driver as 
a tester. Probing and visualisation were carried out us- 
ing the ECL Direct Nucleic Acid Labelling and Detec- 
tion System (Amersham) according to the manufac- 
turer's instructions. Those clones which were positive 
for (a) but negative for (b), or showed a substantially 
larger positive signal with (a) compared to (b), were 
chosen for further analysis. 

DNA sequencing 

Positive clones as identified above were sequenced on 
an automated ABI DNA sequencer (Applied Biosys- 
tems, Warrington, UK). 
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Fi£. / : (i4) Subtractive ddRT-PCR patterns obtained from rat liver following 3-day treatment with WY- 14,643 or phenobarbital. Lane 
1, 1 kb ladder, lane 2, genes up-regulated following Wy, 14-643 treatment; lane 3, genes down-regulated following 
Wy, 14-643 treatment; lane 4, genes up-regulated following phenobarbital treatment; lane 5, genes down-regulated following 
phenobarbital treatment; and lane 6, lkb ladder. (B) Subtractive ddRT-PCR patterns obtained from rat liver showing relative 
changes when phenobarbital treated mRNA is subtracted from Wy-H^i-treated mRNA and vice- versa. Lane 1, 1 kb 
ladder, lane 2, genes showing increased expression following Wy- 14,643 treatment compared to phenobarbital treatment; 
lane 3, genes showing increased expression following phenobarbital treatment compared to Wy- 14,643 treatment. See 
Materials and Methods for further details. 




Fig. 2 : Re-amplified ddRT-PCR products which were down-regulated following phenobarbital treatment (upregulated bands were also 
re-amplified but gel not shown). Individual DNA bands excised from gel of ddRTR-PCR reactions were extracted, 
re-amplified and run on agarose gels to confirm amplification of correct band (numbered). See Materials and Methods for 
further details. 
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Table I : Rat liver genes down-regulated by phenobarbital treatment 

Phenobarbital down-regulated 



Band number (Fig. 2) 
(Approximate size in bp) 



Highest sequence homology 



FASTA-EMBL gene identification 



1 (1500) 




95.3% 


Rat mRNA for 3-oxoacyl-CoA thiolase 


2 (1200) 




92.3% 


Rat hemopoxin mRNA 


3 (1000) 




91.7% 


R rattus alpha-2u-globulin mRNA 


7 (700) 


Clone 1 


77.2% 


M. musculus mRNA for CI inhibitor 




Clone 2 


94.5% 


Rat electron transfer flavoprotein 




Clone 3 


91 .0% 


Mouse topoisomerase 1 (Topo 1) mRNA 


8(650) 


Clone 1 


86.9% 


Scares 2NbMT M. musculus (EST) 




Clone 2 


96.2% 


Rat alpha-2u-globulin (s-type) mRNA 


9(600) 


Clone 1 


86.9% 


Soares mouse NML M. musculus (EST) 




Clone 2 


82.0% 


Soares p3NMF19.5 M. musculus (EST) 


10 (550) 




73.8% 


Soares mouse NML M. musculus (EST) 


11 (525) 




95.7% 


NCLCGAP_Pr1 H. sapiens (EST) 


1*2 (375) 




100.0% 


ft norvegicus mRNA for ribosomal protein 


13 (230) 


Clone 1 


97.2% 


Soares mouse embryo NbME135 (EST) 




Clone 2 


100.0% 


Rat fibrinogen B-beta-chain 




Clone 3 


100.0% 


Rat apoiipoprotein E gene 


14 (170) 




96.0% 


Soares p3NMF19.5 M. musculus (EST) 


15 (140) 




97.3% 


Stratagene mouse testis (EST) 


Others: (300) 




96.7% 


R norvegicus RASP 1 mRNA 


(275) 




93.1% 


Soares mouse mammary gland (EST) 



EST = expressed sequence tag. 
Bands 4-6 were shown to be false positives by dotblot analysis and, therefore, not sequenced. 



Table II : Rat liver genes up-regulated by phenobarbital treatment 



Band number 
(Approximate size in bp) 



Phenobarbital up-regulated 



H ighest sequence homology 



FASTA-EMBL gene identification 



5 (1300) 
7 (1000) 

8(950) 

10 (850) 

11 (800) 

12 (750) 

15 (600) 

16 (550) 
21 (350) 



Clone 1 
Clone 2 



Clone 1 
Clone 2 



93.5% 
95.1% 

98.3% 
95.7% 
94.9% 
75.3% 
93.8% 

92.9% 

95.2% 
93.6% 
99.3% 



Rat cytochrome P450IIB1 

mRNA for rat preproalbumin 

Rat serum albumin mRNA 

NCLCGAP_PM H. sapiens (EST) 

Rat cytochrome P450IIB1 

Rat cytochrome P450IIB1 

Rat cytochrome p450-L (p450IIB2) 

Rat TRPM-2 mRNA 

Rat mRNA for sulfated glycoprotein 

mRNA for rat preproalbumin 

Rat serum albumin mRNA 

Rat cytochrome P450IIB1 

Rat haptoglobulin mRNA partial alpha 

ft norvegicus genes for 18S, 5.8S & 28S rRNA 



EST = expressed sequence tag. 
Bands 1-4. 6, 9, 13, 14 and 17-20 shown to be false positives by dotblot analysis and. therefore, not sequenced. 
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Identification of differentially-regulated 
genes 

* Gene-sequences were identified using the FASTA pro- 
gramme (http^/www.ebi.ac.uk/htbin/fasla.py?request) 
to search all EMBL databases for matching DNA se- 
quences. 

RESULTS 

Figure 1A,B shows the ddRT-PCR patterns of genes 
showing altered expression in rat liver following 3 day 
treatment with phenobarbitaJ or Wy-14,643. Individual 
bands were isolated from the phenobarbitai-modulated 
patterns (both up- and down-regulated), re-amplified 
(Fig. 2), cloned, screened for false positives and then 
identified. Those xenobiotic-modulated gene products 
identified to date are listed in Tables I and II. 

DISCUSSION 

The advent of combinatorial chemistry has led to the 
synthesis of millions of new chemical compounds, 
many of which may be potentially useful in pharma- 
ceutical, agricultural or industrial applications. How- 
ever, whilst there are tests available for those posing a 
genotoxic activity, there remains no short-term assay 
able to identify those chemicals which may belong to 
the non-genotoxic group of carcinogens. 

We have used an adaptation of the subtractive hy- 
bridisation method - ddRT-PCR - to produce charac- 
teristic profiles or 'fingerprints' of those genes which 
are up-regulated or down -regulated in male rat liver 
following acute exposure to test chemicals. The ddRT- 
PCR profiles are characteristic and unique for each of 
the 2 compounds studied to date. 

A number of those gene species showing altered 
expression following phenobarbitaJ treatment have 
been cloned and identified (Tables I & II). It is inter- 
esting to note the presence of CYP2B2 in the up-regu- 
lated genes. This would, of course, be expected fol- 
lowing exposure to phenobarbital and serves as a posi- 
tive control for the method. Other genes which one 
might normally expect to be up-regulated do not ap- 
pear in the table. However, it should be noted that not 
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all bands seen on the differential display were ex- 
tracted and re-amplified due to their being too faint or 
too close to other bands to accurately excise. Further* 
more, it has been well documented [(5) and references 
therein] that a single band extracted from a differential 
display often represents a composite of heterogeneous 
products. We are currently examining new methods to: 
(i) improve resolution of the differential display pat- 
terns (including 2-D agarose gels); and (ii) distinguish 
those ddRT-PCR products which are identical in size, 
but different in sequence. 

Our future efforts will be directed towards deter- 
mining the extent of modulation of a number of the 
genes reported herein using semi-quantitative RT- 
PCR. This should reveal the extent of changes in ex- 
pression of key gene products which may be involved 
in non-genotoxic hepatocarciriogenesis and thus help 
increase understanding of this process. Furthermore, it 
is anticipated that aligning ddRT-PCR profiles of dif- 
ferent non-genotoxic agents found in responsive and 
non-responsive species may enable identification of 
those genes which are mechanistically relevant to the 
non-genotoxic hepatocarcinogenic process. Accord- 
ingly, this approach lends itself well to the identifica- 
tion, characterisation and sub-classification of possible 
different classes of non-genotoxic carcinogens. 
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Abstract 

Understanding the genetic profile of a cell at all stages of normal and carcinogenic development should provide an 
essential aid to developing new strategies for the prevention, early detection, diagnosis and treatment of cancers. We 
have attempted to identify some of the genes that may be involved in peroxisome-proliferator (PP)-induced 
non-genotoxic hepatocarcinogenesis using suppression PCR subtractive hybridisation (SSH). Wistar rats (male) were 
chosen as a representative susceptible species and Duncan-Hartley guinea pigs (male) as a resistant species to the 
hepatocarcinogenic effects of the PP, [4-chloro-6-(2,3-xylidino)-2-pyrimidinylthio] acetic acid (Wy-1 4,643). In each 
case, groups of four test animals were administered a single dose of Wy-1 4,643 (250 mg/kg per day in corn oil) by 
gastric intubation for 3 consecutive days. The control animals received corn oil only. On the fourth day the animals 
were killed and liver mRNA extracted. SSH was carried out using mRNA extracted from the rat and guinea pig 
livers, and used to isolate genes that were up and downregulated following Wy-1 4,643 treatment. These genes 
included some predictable (and hence positive control) species such as CYP4A1 and CYP2C11 (upregulated and 
downregulated in rat liver, respectively). Several genes that may be implicated in hepatocarcinogenesis have also been 
identified, as have some unidentified species. This work thus provides a starting point for developing a molecular 
profile of the early effects of a non-genotoxic carcinogen in sensitive and resistant species that could ultimately lead 
to a short-term assay for this type of toxicity. © 2000 Elsevier Science Ireland Ltd. All rights reserved. 

Keywords: Wy-1 4,643; Peroxisome proliferator; Non-genotoxic hepatocarcinogenesis; Suppression PCR subtractive hybridisation; 
RT-PCR; Rat; Guinea pig; Gene regulation; Differential gene display; Gene profiling 
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Introduction 

rhe advent of combinatorial chemistry and 
nputer-aided drug design has led to a recent 
surge in the number of chemical compounds 
it have potential therapeutic, agricultural and 
lustrial applications. Although it has been sug- 
ted that the contribution of synthetic chemicals 
the overall incidence of human cancer is low, 
re still remains an absolute requirement to 
tluate all new chemicals for toxic and carcino- 
id potential. The latter is one of the most 
>blematic areas of chemical safety evaluation 
1 is usually carried out using short-term in vitro 
1 in vivo genotoxicity assays augmented by 
onic bioassay tests. The short-term assays have 
tved useful in the early identification of poten- 
genotoxic carcinogens, but their value is lim- 
I by observations that suggest that 
>roximately 60% of chemicals identified as car- 
;>gens in life-exposure studies produce mainly 
ative findings in short-term genotoxcity tests 
hby, 1992; Parodi, 1992). Thus, there is cur- 
tly no reliable and rapid means of evaluating 
carcinogenic risk of new .chemicals that fall 
> this latter group of compounds, termed non- 
otoxic (or epigenetic) carcinogens. 
)ne approach to addressing this problem is to 
;idate the molecular mechanisms by which 
<wn non-genotoxic carcinogens act. It should 
l be possible to identify common factors/ 
:hanisms that can serve as early biomarkers of 
;inogenic potential for new chemicals. To this 
, a large number of groups have reported on 
various effects of non-genotoxic compounds 
various animal species (Marsman et al., 1988; 
e et al., 1993; Cattley et al., 1994; Hayashi et 
1994; Human and Experimental Toxicology, 
1; Anderson et al., 1996). However, the mech- 
tic picture is still far from complete with many 
hose genes involved in the carcinogenic pro- 
remaining unknown, and their identification 
efore remains a key goal in elucidating the 
ecular mechanisms by which non-genotoxic 
inogenesis occurs. 

jbtractive hybridisation (SH) and related tech- 
>gies such as representational difference analy- 
(RDA) (Hubank and Schatz, 1994) and 



differential display (DD) (Liang and Pardee 
1992) can be used to aid the isolation of genes 
showing altered expression in target tissues fol- 
lowing exposure to a chemical stimulus. These 
techniques can also be used to identify differential 
gene expression in neoplastic and normal cells 
(Liang et al., 1992), infected and normal cells 
(Duguid and Dinauer, 1990), differentiated and 
undifferentiated cells (Sargent and Dawid, 1983; 
Guimaraes et al., 1995), activated and dormant 
cells (Gurskaya et al., 1996; Wan et al., 1996), 
different cell types (Hedrick et al., 1984; Davis et 
al., 1984) amongst others. Most importantly, us- 
ing such approaches, no prior knowledge of the 
specific genes that are upregulated/downregulated 
is required. 

Using a variation of SH, termed suppression- 
PCR subtractive hybridisation (SSH) (Diatchenko 
et al., 1996), we have previously reported the 
isolation of a number of genes showing altered 
expression in male rat liver following acute expo- 
sure to phenobarbital (Rockett et al., 1997). In 
the current work we have used the same experi- 
mental approach to isolate genes that are differen- 
tially expressed in the livers of male rats and 
guinea pigs following short-term (3-day) exposure 
to the peroxisome proliferator (PP) and non- 
genotoxic hepatocarcinogen, Wy-14,643. We have 
isolated and identified a number of gene species, 
some of which may be important in the induction 
of, or protection against, non-genotoxic 
hepatocarcinogenesis. 



2. Materials and methods 

2.7. Animals and treatment 

All animal experiments were undertaken in ac- 
cordance with Her Majesty's Home Office De- 
partment guidelines under the auspices of 
approved personal and project licences. Male 
Wistar rats (150-200 g) and male Duncan-Hart- 
ley guinea pigs (250-300 g) were obtained from 
Kingman and Bantam (Hull, UK). Upon receipt, 
both groups were randomly assigned into two 
groups of four. They were maintained on a rat, 
mouse or guinea pig standard diet (B&K Univer- 
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sal, Hull) and a daily cycle of alternating 12-h 
. periods of dark and light. The room temperature 
was maintained at 19°C and a relative humidity of 
55%. The animals were acclimatised to this envi- 
. ronment for 7 days before treatment commenced. 
[4-chloro-6-(2,3^ylidino)-2-pyrimidinylthio] acetic 
acid (Wy-14,643, Campo, Emmerich; 250 mg/kg 
per day in corn oil) was administered by gavage 
to the treated groups of rats and guinea pigs on 3 
consecutive days, whilst control groups received 
an equal volume of corn oil only. During this 
time, all animals had free access to food and 
water. The animals were killed by cervical disloca- 
tion on the fourth day, and their livers immedi- 
ately excised, weighed, sliced into approximately 
0.5-cm cubes, snap frozen in liquid nitrogen and 
stored at - 70°C. 

2.2. mRNA extraction 

Approximately 0.25 g of each frozen liver sam- 
ple was ground under liquid nitrogen using a 
mortar and pestle. Messenger RNA was extracted 
from the ground liver using the PolyATtract® 
System 1000 kit (Promega, Madison, USA) ac- 
cording to the technical manual provided by the 
manufacturers. The mRNA was DNase-treated 
(RQ Rnase-free Dnase, Promega, final concentra- 
tion 10 U/ml) before phenol/chloroform extrac- 
tion and ethanol precipitation. The mRNA was 
redissolved at a final concentration 500-1000 ng/ 
Ml 

23. cDNA Subtraction 

This was carried out using the PCR-Select™ 
cDNA Subtraction Kit (Clontech, Palo Alto, 
USA) according to the manufacturer's instruc- 
tions. Subtractions were carried out with mRNAs 
derived from single animals. The mRNA from the 
remaining three animals in each group was later 
used for quantitative RT-PCR analysis of specific 
genes. 

2.4. Band extraction and cloning 

The secondary PCR reactions from the cDNA 
subtraction procedure were run on a 2% 
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Metaphor agarose gel (FMC, Rockland, USA) 
containing 0.5 ng/ml ethidium bromide (Sigma, 
Dorset, UK). One times TAE (0.04 M Tris-ac- 
etate, 0.001 M EDTA) was used to prepare the gel 
and as the running buffer. After running for 6-7 
h at 3.75 V/cm, the gel was overstated for 30 min 
with SYBR Green I DNA stain (FMC, 1:10000 
dilution in 1 x TAE). Each discernible band of 
the differential display pattern was extracted from 
the gel with a scalpel and the DNA eluted using a 
Genelute™ agarose spin column (Supelco, Belle- 
fonte, USA). Five microlitres of the eluted DNA 
was reamplified using the original nested (sec- 
ondary) PCR primers supplied with the PCR-Se- 
lect™ cDNA subtraction kit. The PCR products 
were electrophoresed on a 2% standard agarose 
gel (Boehringer Mannheim, East Sussex, UK) and 
the reamplified target bands extracted from the 
gel as above. The eluted DNA was immediately 
ligated into a TOPO TA Cloning® vector (Invitro- 
gen, Carlsbad, USA) before transformation in 
Escherichia coli TOPI OF 9 One Shot™ cells 
(Invitrogen). 

2.5. Colony screening 

2.5 A. Stage I 

Eight transformed (white) colonies from each 
band were grown up for 6 h in 200 |il LB broth 
containing ampicillin (Sigma, 50 mg/ml). One mi- 
crolitre of this was subjected to PCR using the 
same conditions and nested primers as described 
above. One tenth (2 \x\) of the completed PCR 
reaction was electrophoresed on a 2% standard 
agarose gel and one tenth on a 2% standard 
agarose gel containing HA red (Hanse Analytik 
GmbH, Bremen, Germany, 1 U/ml) to discern the 
differentially cloned products. The remainder of 
the PCR reaction was used to prepare duplicate 
dotblots on Hybond N + membranes (Amersham, 
Little Chalfont, UK). 

2.5:2. Stage II 

The duplicate dotblots were probed with (a) the 
final differential display reaction and (b) the 're- 
verse-subtracted' differential display reaction. To 
make the 'reverse-subtracted* probe, the subtrac- 
tive hybridisation step of the differential display 
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T-PCR procedure was carried out using the 
riginal tester (treated) mRNA as the driver and 
{e original driver (control) mRNA as the tester, 
robing and visualisation were carried out using 
le ECL direct nucleic acid labelling and detec- 
on system (Amersham, Little Chalfont, UK) ac- 
xding to the manufacturer's instructions. Those 
ones that were positive for (a) but negative for 
)), or showed a substantially larger positive sig- 
al with (a) compared to (b), were selected for 
>NA sequence analysis. 

6. DNA sequencing 

The remainder of the cultures (prepared in 
age 1 screening) containing different cloning 
roducts (as discerned in the two screening steps) 
ere grown up overnight in 5 ml LB broth con- 
.ining ampicillin (50 mg/ml). A plasmid miniprep 
as prepared from each (Wizard Plus SV 
[inipreps DNA purification system, Promega) 
xording to the manufacturer's instructions. The 
oned inserts were sequenced on an automated 
BI DNA sequencer (Applied Biosystems, War- 
ngton, UK) using the M13 forward primer 
jTAAAACGACGGCCAGT) or M13 reverse 
imer (AACAGCTATGACCATG). 

7. Identification of differentially regulated genes 

Gene sequences thus obtained were identified 
;ing the FASTA 3.0 programme (Lipman and 
jarson, 1985; Pearson and Lipman, 1988) (http:/ 
/ww.ddbj.nig.ac.jp/E-mail/homology.html) to 
arch all EMBL databases for matching DNA 
quences. Each clone sequence was submitted in 
e forward and reverse direction, and the one 
turning the highest statistical probability of 
atch to a known sequence was noted. Sequence 
jmologies between our submitted clone sequence 
id the queried database sequence were deter- 
ined (by FASTA) over a region of at least 60 
ise pairs. 

S. RT-PCR analysis of selected candidate genes 

cDNA sequences of the target genes were ob- 
med from the NIH gene database (GenBank at 
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http: / / www.ncbi.nlm.nih.gov/Web / Search / index, 
html) and the computer programme gene 
jockey (BioSoft, Cambridge, UK) used to select 
primer pairs from these sequences. Where guinea 
pig sequences were available, rat and guinea pig 
sequences were aligned and primers chosen from 
regions of homology. If guinea pig sequences were 
not available, rat and human sequences were 
used. In cases where exact homology could not be 
found, the sequence from the rat was used. In the 
case of CD81 only, no rat or guinea pig sequences 
were available and so mouse and human se- 
quences were aligned and a primer pair chosen 
from a region of homology. Primers (obtained 
from Gibco-BRL, Paisley, UK) were dissolved at 
a concentration of 50 pmol/jal in sterile distilled 
water and stored at — 20°C The primer pairs 
used plus other reaction parameters are shown in 
Table 1. mRNA was extracted (as described 
above) from all four treated animals and from 
three animals in the control group. Integrity of 
the eluted mRNA was confirmed on a 2% agarose 
gel, and the concentration and purity were mea- 
sured using a Genequant II spectrophotometer 
(LKB, Bromma, Sweden) and then diluted to 10 
ng/jil. One microlitre of this latter solution was 
used per RT-PCR reaction. 

RT-PCR was carried out in a single tube (50 nl) 
reaction using the Access RT-PCR system 
(Promega) according to manufacturer's instruc- 
tions. In the kinetic and quantitative analyses, 
omission of RNA was used as a control for the 
presence of any contaminating DNA. After ob- 
taining a PCR signal of the correct size and 
optimising the reaction conditions, each PCR 
product was digested with between two and four 
separate restriction enzymes. Specific restriction 
patterns were thus obtained, which further confi- 
rmed the identity of the PCR products as being 
the original target genes. Kinetic analysis (14-32 
cycles) was then performed in each case to deter- 
mine the location of the mid-log phase. 

For the semi-quantitative analysis of each 
target gene, RT-PCR reactions were carried out in 
triplicate for each sample to reduce the effect of 
intertube RT-reaction variations (Kolls et al., 
1993) and pipetting errors. For each gene, a mas- 
termix containing enough reagents for three times 
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the number of samples (seven for rat, six for 
guinea pig) was prepared except that mRNA was 
omitted, the latter being added after aliquoting 49 
^1 of the mastermix into an appropriate number 
of tubes. Amplification of albumin (the reference 
gene) was carried out in separate tubes since the 
mid-log phase of this gene is at a much lower 
cycle number than the target genes due to its high 
abundance. All RT-PCR products were analysed 
on 2% agarose gels containing 0.5 |ig/ml ethidium 
bromide. The target gene samples were loaded on 
the gel first and run in at 3 V/cm for 10 min. The 
corresponding albumin samples were then loaded 
and the gel run for a further 1/2 h. In this way, all 



RT-PCR products from each target gene and 
albumin from the corresponding samples could be 
run on the same gel. Gels were photographed 
using type 665 posi-neg film (Sigma) and quanti- 
tation of the band intensity was carried out using 
a dual wavelength flying spot laser scanner densit- 
ometer (Shimadzu). 

2.9. Statistical analysis 

Statistical analysis of unpaired samples was car- 
ried out using the two-tailed Student's f-test. Val- 
ues were considered statistically significant at 
P < 0.05 or less. 



L 1 2 



L 1 2 





-. A-.'.".. 



m 

• Ay, ■,;>, 



B 

Fig. I. Final displays of differentially expressed genes that 
were (I) upregulated and (2) dowhregulated in rat (A) and 
guinea pig (B) livers following 3-day treatment with Wy- 
14,643. mRNA extracted from control and treated livers was 
used to generate the differential displays using the PCR-Select 
:DNA subtraction kit (Clontech). Lane (L) is a 1 Kb DNA 
Ladder standard and 10 ul of secondary PCR reaction were 
loaded in all other lanes. 



3. Results 

3.1. Cloning and screening of transcripts 

For both the rat and guinea pig experimental 
groups, cDNA subtraction was carried out in the 
forward (control driving tester) and reverse (tester 
driving control) directions to isolate both upregu- 
lated and downregulated mRNA species respec- 
tively. Using a standard primary hybridisation 
time of 8 h we obtained a substantial amount of 
non-specific products in all the final differential 
displays (data not shown). This background 
smearing was almost completely removed by re- 
ducing the primary hybridisation time to 4 h 
(CLONTECHniques, 1996). Fig. 1 shows the 
ddRT-PCR patterns of genes showing altered ex- 
pression in rat and guinea pig liver following 
3-day treatment with Wy- 14,643. The profiles are 
unique for each species, and in each case the 
profile for the upregulated genes (control mRNA 
driving tester mRNA) is different to that obtained 
for the downregulated genes (tester mRNA driv- 
ing control mRNA). 

The practical outcome of the SSH method is 
that a series of differentially expressed genes is 
observed as a ladder on an agarose gel. The 
majority of these gene fragments fall within the 
150-2000 bp range, with bands up to 5 Kbp 
occasionally being observed. Each band may the- 
oretically consist of one or more products of 
similar size, as the gel has a maximum resolution 
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Fig. 2. Discrimination of different ddRT-PCR products having 
the same molecular size using HA-red. Gel (A) is a 2% 
standard agarose gel Gel (B) is a 2% standard agarose gel 
containing 1 U/ml HA-red. Band numbers refer to the sequen- 
tial bands (largest to smallest) extracted from the original 
display of genes upregulated in rat liver following 3-day treat- 
ment with Wy-14,643. Ten micorlitres of each PCR reaction 
were loaded per lane. 



of approximately 1.5% (3 bp per 200). In addi- 
tion, there may be two or more products that are 
the 'same size, but have a different sequence. 



Therefore some form of discrimination must be 
employed to isolate as many of these products as 
possible. HA-red screening (Geisinger et ah, 1997) 
of a number of clones derived from each band 
provided a means to discriminate between differ- 
ent gene species of the same size. A typical exam- 
ple of such a gel is shown in Fig. 2. In total, 88 
and 48 apparently different clones were obtained 
from the final differential expression patterns of 
upregulated and downregulated rat genes, respec- 
tively. Sixty nine and 89 apparently different 
clones were obtained from the final differential 
expression patterns of the upregulated and down- 
regulated guinea pig genes, respectively. 

Having identified as many different candidate 
gene products as possible in the screening step I, a 
second screening step was carried out on every 
clone to confirm those that represented true dif- 
ferentially expressed genes. This is necessary since 
no subtraction technique is 100% efficient. The 
approach we used, termed PCR-select differential 
screening (as recommended in Clontech's PCR-se- 
lect cDNA subtraction, kit protocol), utilises the 
forward and reverse subtractions as an aid to 
screening for the true differentially expressed 
genes (CLONTECHniques, 1997). Because these 
probes have already undergone subtraction, they 
have been enriched for differentially expressed 
genes and are therefore more sensitive than un- 
subtracted driver/tester cDNA probes for detect- 
ing true differential expression. All the clones that 
were isolated from each display were dotblotted 
and probed with the display from which they was 
obtained, plus the corresponding reverse-sub- 
tracted display. An example of such a blot is 
shown in Fig. 3. Clones corresponding to authen- 
tic differentially expressed mRNAs hybridised 
with the subtracted cDNA probe, but not the 
reverse-subtracted probe. We also included in the 
authentic positives, those clones that gave a sub- 
stantially greater signal with the subtracted probe 
compared to the reverse-subtracted probe. False 
positives hybridised with either both probes or 
with neither probe. Of the original 88 upregulated 
and 48 downregulated rat clones selected for this 
screening step, 28 (32%) and 15 (31%) respec- 
tively, were found to be true positives. In the rat, 
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! (100%) of the true positive upregulated genes 
able 2) and 1 1 (73%) of the true positive down- 
gulated genes (Table 3) were non-redundant. Of 
e original 69 upregulated and 89 downregulated 
linea pig clones selected for this screening step, 
! (70%) and 37 (42%) respectively, were found to 
; true positives. Thirty six (75%) of the upregu- 
ted genes (Table 4) and 33 (89%) of the down- 
gulated genes (Table 5) were non-redundant. 

2. Identification of clones 

On sequence analysis it was found that some 
imes were unsequencable in the first instance 
113 forward primer) due to long polyA runs 
at appeared to prematurely terminate the se- 
lencing reaction. These clones were therefore 
sequenced from the opposite direction using the 
13 reverse primer. Those xenobiotic-modulated 
ne products identified to date are listed in Ta- 
*s 2 and 3 (rat) and Tables 4 and 5 (guinea pig). 




. 3. Dot blots of clones of putative upregulated gene species 
ated from guinea pig liver following 3-day treatment with 
-14,643. AH clones identified in the stage I screening step 
methods) were blotted and probed with (A) the differen- 
display from which they originated (control driving 
ted) and (B) the reverse subtraction (treated driving con- 
). Arrows indicate some of the true differentially expressed 
les. 



Table 2 

Identification of genes that were upregulated in male rat liver 
following 3-day treatment with WY- 14,643 



FASTA-EMBL gene Accession No. Sequence 
identification (rat un- homology* (%) 

less otherwise stated) 
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Table 2 (Continued) 



FASTA-EMBL gene 
identification (rat un- 
less otherwise stated) 



Accession No. 



Sequence 
homology 0 (%) 



17-P-hydroxsteroid de- 
hydrogenase 

Soares mouse 
p3NMF19.5 (EST) 

Peroxisomal enoyl- 
CoA:hydratase -3- 
hydroxyacyl CoA 
bifunctional enzyme 

Integral membrane 
protein, TAPA-1 
(CD81) (mouse) 

Soares mouse lymph 
node (EST) 

H. sapiens (clone 
zapl28) mRNA 

Lysophospholipase ho- 
mologue (human) 

Soares mouse lymph 
node (EST) 



RN17BHDT2 

AAO38051 

RNPECOA 



87 
87 
85 



S45012 

MMAA88445 
L40401 
HSU67963 
AA2 1 7044 



81 

81 
76 
76 
74 



0 Refers to the nucleotide sequence homology between the 
cloned band isolated from the differential display and the 
corresponding gene derived from the EMBL gene sequence 
bank. 

b EST is 'expressed sequence tag' — a gene of as yet 
unknown identity and function. 

c Where sequence homologies were equal in both directions 
of the isolated band, both the sense (S) and antisense (A) 
identities are given. 



Table 3 

Identification of genes that were downregulated in male rat 
liver following 3-day treatment with Wy-14,643 



FAST-EMBL gene Accession No. Sequence 
identification (rat un- homology" (%) 

less otherwise stated) 



NCI_CGAP_Lil (H. 
sapiens) (EST b )(S c ) 

NCI_CGAP_Prl (H. 
sapiens) (EST)(AS C ) 

UDP-glucuronosyl- 
transferase 
(UGT2B12) 

Complement compo- 
nent c3 

Soares mouse pla- 
centa (S) 

Ape (chimpanzee) 28S 
rRNA (AS) 

Rat CYP2C11 

Ribosomal protein S5 

Transthyretin 

Contrapsin-like 
protease inhibitor 

Prostaglandin F2a (S) 

P-2-microglobulin 
(AS) 

Apolipoprotein C-III 
Parathymosin-alpha 

(zinc2 + -binding 

protein) 



AA484528 99 

AA469320 99 

RN06273 98 

RNC3 96 

AA023305 96 

PTRGMC 96 

RNCYPM 1 95 

RNRPS5 94 

RNTTHY 94 

RNCCP23 89 

RN26663 84 

RNB2MR 84 

RNAPOA02 82 

RN11ZNBP 75 



" Refers to the nucleotide sequence homology between the 
cloned band isolated from the differential display and the 
corresponding gene derived from the EMBL gene sequence 
bank. 

b EST is 'expressed sequence tag* — a gene of as yet 
unknown identity and function. 

c Where sequence homologies were equal in both directions, 
both the sense (S) and antisense (A) identities are given. 



In all cases, both the forward and reverse se- 
quence of the target clones were analysed and the 
gene having the highest statistical homology 
noted; 

33. RT-PCR analysis of selected clones 

The results of a typical RT-PCR semi-quantita- 
tion experiment for transferrin in the rat is given 
in Fig. 4 and the results for a total of 12 selected 
genes in both the rat and guinea pig are shown in 
Table 6. 



4. Discussion 

It is now apparent that all cancers arise from 
accumulated genetic changes within the cell. Al- 
though documenting and explaining these changes 
presents a formidable obstacle to understanding 
the different mechanisms of carcinogenesis, the 
experimental methodology is now available to 
begin attempting this difficult challenge. In order 
to beein the elucidation of the molecular rnecha- 
hisms involved in non-genotoxic hepatocarcino- 
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lesis, we have used SSH to identify a number of 
les that are upregulated or downregulated in 
tie rat and guinea pig livers following short 
m exposure to the PP, Wy- 14,643. We have 
id the rat model to represent a species suscepti- 
; to the non-genotoxic carcinogenic effect of 
•s and the guinea pig as a resistant species 
rton et al., 1984; Rodricks and Turnbull, 1987; 
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Lake et al., 1989; Makowska et al., 1992; Lake et 
al, 1993). 

Gurskaya et al. (1996), who originally devel- 
oped the SSH technique, cloned the products of 
the secondary PCR reaction and screened a small 
number of randomly selected colonies for differ- 
entially expressed clones using northern hybridisa- 
tion. However, we decided against this approach 



>Ie 4 

atification of genes that were upregulated in male guinea pig liver following 3-day treatment with WY-14,643 

STA-EMBL gene identification (guinea pig unless otherwise stated) Accession No. Sequence 

homology" (%) 
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AA 156847 
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80 


oplasmic chaperonin hTRiCS (human) 


U17104 


80 


ha-l-antiproteinase F 


M57271 


77 


erogeneous nuclear ribonuclearprotein cl/c2 (human) 


D28382 


77 


res parathyroid tumour (EST) (similar to human serum albumin precursor) 


AA860651 


76 


itagene mouse kidney (EST) 


AA 107327 


75 


res parathyroid tumour NbHPA human cDNA (EST) 


AA860653 


74 


res mouse mammary gland (EST) 


AA6 19297 


74 


NA clone 15 004 (EST) (human) 


H01826 


74 


res senescent fibroblasts (EST) (mouse) 


W52190 


74 


proalbumin (human) 


E04315 


72 


NA clone 73 169 (EST) (human) 


T56624 


72 


imin D-binding protein (human) 


L10641 


71 


>H gene (exon 8) (human) 


Y 11498 


71 


IL flow sorted chromosome 


B05457 


71 


res foetal liver spleen (EST) (mouse) 


AA009524 


71 


res foetal heart NbMH19W (EST) (mouse) 


AA009421 


69 


res foetal heart NbHH19W H. sapiens cDNA clone (EST) 


W94377 


67 


nylalanine hydroxylase (human) 


U49897 


67 


line-5-carboxylate dehydrogenase (human) 


U24266 


66 


tathione-5-transferase homologue (human) 


U90313 


65 


l_CGAP_GCBI (EST) (human) 


AA769294 


65 


* 

tective protein (human) 


M22960 


64 


nc 27 375 (EST) (human) 


N37046 


62 


tasene colon ( # 937 204) H, sapiens cDNA clone (EST) 


AA149777 


62 



Refers to the nucleotide sequence homology between the cloned band isolated from the differential display and the correspond- 
gene derived from the EMBL gene sequence bank. 
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.Table 5 

. Identification of genes that were downregulated in male guinea 
pig liver following 3-day treatment with WY-14,643 

FASTA-EMBL gene Accession No. Sequence 
- identification (guinea homology* (%) 

pig unless otherwise 
stated) 



Complement 


M34054 


97 


• 

protein 




95 


Munnoglobuhn 


Alpha- 1 -an- 


MS7271 


88 


tiproteinase F 


AO**"-' 


89 


Elongation tactor-ai- 


pha-1 (rabbit) 






Coupling protein O 




oo 


(human) 






NCI_CGAP_Ovl 




R7 
o / 


(EST b ) (human) 




O J 


Lecithin xholesterol 


acetyl transierase 






(rabbit) 




84 


All 1 Tl \ 

Aldolase B (human) 


Anti-thrombin III 


E00116 


80 


(human) 






Phenylalanine ny- 






droxylase (human) 






Inter-a-trypsin in- 




70 


hibitor (human) 






Normalised rat mus- 




7R 


cle (Jbo l ) (o ; 






Normalised rat ovary 




"75 








Complement factor 




11 


Ba fragment (hu- 






man) 






Dihydrodiol denydro- 




76 


genase (human) 


Y08409 


75 


opotl** gene ^inyruio- 


inducible hepatic 






proteinxmJtnanj 






HAL. clone i /**pix 


AP004236 


75 


(human) 






Mitochondrial alde- 


X05409 


74 


hyde dehydroge- 






nase (human) 






Preproalbumin (hu- 


E04315 


74 


man) 






NCI_CGAP_Pr9 


AA533142 


74 


(EST) (human) (S) 






Normalised rat pla- 


AA851197 


74 


centa (EST) (AS) 






Heparin sulfate pro- 


J04621 


73 


teoglycan (human) 






cDNA clone 33 992 


R24330 


73 


(EST) (human) 







Table 5 (Continued) 



FASTA-EMBL gene Accession No. 


Sequence 


identification (guinea 


homology* (%) 


pie unless otherwise 




stated) 




Retinoi dehydrogenase U3350I 


71 


(rat) 




TAPA-1 integral mem- S45012 


71 


hrane nrotein 




(CT*>8n (mouse! 






70 


JJClll 




Annlinonrntein B (t>\£\ LI 1235 


69 


rDNA clone 143 918 R76742 


68 






rr-fiHrinr>pen (human! K.02569 


68 


Qr»ar#»c foptnl livpr Vym736 


68 


^nleen INF (mouse! 




Raroead bowel (EST! AA232049 


67 








66 


iransi erase icaij 




IViyeiOJu ICUJtaCnila CCIl LUOZHU 


65 


aiiiercriijdijon 








(human) (S) 




STS SHGC-34 987 (hu-G27984 


65 


man) (AS) 




Soares mouse AA222798 


64 


3NME125 




Strata gene mouse em- AA 199420 


64 


bryonic (EST) (S) 




Rad 52 (mouse) AF004854 


63 



a Refers to the nucleotide sequence homology between the 
cloned band isolated from the differential display and the 
corresponding gene derived from the EMBL gene sequence 
bank. 

b EST is 'expressed sequence tag' — a gene of as yet 
unknown identity and function 

c Where sequence homologies were equal in both directions, 
boththe sense (S) and antisense (A) identities are given. 

for several reasons: (1) the kinetics ofligation and 
transformation favour the isolation of smaller 
PCR products, thereby producing a misrepresen- 
tation of larger gene products; (2) northern blot 
analysis is notoriously insensitive and is unlikely 
to confirm expression of rare transcripts; (3) there 
is no measurable end point to the screening of 
clones produced in this way other than to analyse 
every transformed colony. We used instead an 
alternative approach; after running out the differ- 
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:ntial display on a high-resolution agarose gel 
Fig. 1) and overstating with SYBR Green I to 
:iihance visualisation, the composite bands were 
ndividually extracted, reamplified and cloned, 
iowever, it has been well documented that single 
>ands from differential displays often contain a 
leterogeneous mixture of different products 
Mathieu-Daude et al'., 1996; Smith et aL, 1997). 
fhis is because polyacrylamide gels cannot dis- 
riminate between DNA sequences that differ in 
ize by less than about 0.2% (Sambrook et aL, 
989). High-resolution agarose gels such as those 
ised in this work are even less sensitive, normally 
>nly discriminating products that differ in size by 
t least 1.5%. The use of the HA-red screening 
tep enables resolution of identical or nearly iden- 
ical sequences based on their AT content (Wawer 
t aL, 1995) and is sensitive down to < 1% differ- 
nce. Furthermore, it is rapid, technically simple 
nd does not require the use of radiolabels. 
}eisinger et al. (1997) originally demonstrated the 
isefulness of using HA-red to identify different 
products cloned from the same band of an RNA 
lifferential display experiment by simultaneously 
unning them in normal agarose (to discriminate 
»y size) and in normal agarose containing HA-red 
to discriminate by AT content). We have found 
hat this approach is equally useful for identifying 
ifferent gene species cloned from the same band 
■f our SSH display. 

Diatchenko et al. (1996) reported that SSH is 
tighly efficient at producing differentially ex- 
pressed gene species. However, we also included a 
econd screening step to further confirm that the 
lones isolated from the differential display were 
adeed differentially expressed. Duplicate dotblots 
■f the candidate clones were blotted with the 
iisplay from which they were originally isolated 
nd with the 'reverse subtraction' display. To 
nake the reverse-subtracted probe, the subtractive 
ybridisation step of the procedure was carried 
■ut using the original tester cDNA as a driver, 
nd the original driver cDNA as a tester. In this 
/ay, clones that are false positives can be iden- 
ified through their presence in both blots. Such 
alse positives most commonly arise through hav- 
ig a very high abundance in the initial sample or 
nusual hybridisation properties (Li et aL, 1994). 



Although the SSH method itself has been 
shown to be efficient, and despite the screening 
step that we included, there is an important caveat 
to bear in mind — namely that it is important 
that all clones be considered only as 'candidates' 
until the actual abundance of their mRNA is 
quantitated in treated and control samples. To- 
wards this end, we examined the expression of a 
limited number of clones using semi-quantitative 
RT-PCR. Albumin was used as the reference eene 
as we have previously found that the expression 
of this gene does not appear to change with the 
treatment regime that we used (Fig. 4, and data 
not shown). There are a number of interesting 
points to note from our results. The first is the 
presence of genes that serve as appropriate posi- 
tive controls in the upregulated and downregu- 
lated series. For example, in the rat it can be seen 
that CYP4AI expression increases 14-fold follow- 
ing treatment. Although CYP4AI mRNA expres- 
sion levels following WY- 14,643 treatment have 
not been previously reported in this model, the 
figure compares favourably with that recorded by 
Bell et al. (1991), who used RNAse-protection to 
quantitate CYP4A1 in rat liver following treat- 
ment with methylclofenapate, another PP. In ad- 
dition, we also confirmed that the peroxisomal 
enoyl-CoA:hydratase-3-hydroxyacyl-CoA Afunc- 
tional enzyme is also upregulated 9-fold, in agree- 
ment with the findings of Chen and Crane (1992). 

A number of genes were downregulated follow- 
ing Wy-14,643 exposure, including CYP2C11 ex- 
pression. Corton et al. (1997) reported similar 
findings and suggested that this may in part ex- 
plain why male rats exposed to Wy-14,643 and 
some other PPs have high serum estradiol levels, 
as estradiol is a substrate for CYP2C11. We have 
also shown that the expression of contrapsin-like 
protease inhibitor (CLPI) was downregulated by 
Wy-14,643. This has not previously been reported, 
and we suggest that it may be linked to a require- 
ment for increased availability of amino acids to 
accommodate the hepatomegaly induced by treat- 
ment. Although little is known of the function of 
parathymosin-a, (zinc 2 + -binding protein) it has 
been shown to interact with the globular domain 
of histone HI, suggesting a role in histone func- 
tion (Kondili et aL, 1996). In contrast to the 
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Albumin 
Transferrin 



Albumin 
Transferrin 




Fig. 4. Semi-quantitative RT-PCR experiment showing relative decrease in expression of transferrin in treated rat liver (RT-1 to 
RT-4) compared to controls (RC-1 to RC-3). An equal amount of mRNA was used in each reaction (10 ng), and each sample was 
quantitated in triplicate to reduce the effects of inter-tube variation. N is negative control (no mRNA). Lane M is a 100 bp ladder 
and lane L is a 1 Kb DNA ladder. 



downregulation observed in this work, other stud- 
ies have shown that parathymosin-ot expression is 
elevated in breast cancer (Tsitsilonis et al., 1993, 
1998), with the implication that parathymosin-a 
may somehow be involved in regulating cell pro- 
liferation by more than one mechanism. Transfer- 
rin has previously been shown to be 
downregulated in rat liver by hypolipidemic PPs 
(Hertz et al., 1996). It is therefore interesting to 
note that we isolated a clone identified as transfer- 
rin from the upregulated display profile. Since we 
confirmed by RT-PCR that transferrin is in fact 
downreeulated in the rat (Fig. 4), we conclude 
that transferrin was either a false positive or was 
incorrectly identified. It could also be that we 
have isolated a close relative, splice variant or 
isoform of transferrin, which demonstrates a dif- 
ferent expression profile under these experimental 
conditions. Further investigations are therefore 



required to determine which of these possibilities 
are correct. 

One of our most intriguing observations was 
that one gene, CD81, appeared to be upregulated 
in rat liver but downregulated in guinea pig liver 
following Wy- 14,643 exposure. CD81 is a widely 
expressed cell surface protein that is involved in a 
large number of cellular functions, including ad- 
hesion, activation, proliferation and differentia- 
tion (reviewed by Levy et al., 1998). Since all of 
these functions are altered to some extent in car- 
cinogenesis, it is perhaps an important observa- 
tion that CD81 expression is differentially 
regulated in a resistant and sensitive species ex- 
posed to a non-genotoxic carcinogen. 

Albumin and ribosomal genes appear common 
to all differential displays and are thus undesir- 
able false positives. However, due to their high 
expression in the liver, they are difficult to re- 
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iove. We also noted a number of gene species, 
articularly in the guinea pig, which were com- 
ion to both upregulated and downregulated 
rofiles. Again, the most likely reason for these 
aving arisen is their high abundance. 
A relatively large number of upregulated and 
ownregulated genes were isolated from guinea 
ig liver following Wy- 14,643 exposure. However, 
le guinea pig genome has been relatively poorly 
haracterised and so many of the clones were 
lentified as resembling genes or ESTs from other 
secies. Without full-length sequence data it is 
ifficult to ascertain the accuracy of the assigned 
lentities and this must be borne in mind when 
tilising data such as this, for example, in design- 
lg effective primers for RT-PCR studies. Al- 
lough the actual isolated clone sequences can be 
sed to do this, their relatively small size often 
jstricts the ability to design effective primers. In 
ddition, as we observed with transferrin, using a 
ublished full-length sequence may help to iden- 
fy false positives. 

able 6 

imi-quantitative RT-PCR analysis of selected gene species in the rat and guinea pig 3 



By comparing the expression profiles of genes 
showing altered expression in a PP-sensitive spe- 
cies (rat) with a PP-resistant species (guinea pig), 
it was our aim to identify genes that are mecha- 
nistically relevant to the non-genotoxic hepatocar- 
cinogenic action of Wy- 14,643. However, few of 
the genes that we have isolated were common to 
both the rat and the guinea pig. This suggests 
either that the molecular mechanisms of response 
in these two species are so different that few genes 
are commonly regulated in response to Wy- 14,643 
exposure, or that we have recovered only a small 
proportion of those genes that have altered ex- 
pression. The latter seems the more likely scenario 
since it is perceived that one of the main problems 
of subtractive hybridisation and other differential 
expression technologies is the inability to consis- 
tently isolate rare gene transcripts (Bertioli et al., 
1995). This is potentially problematic in that 
weakly expressed genes may play an important 
role in regulating key cellular processes, and that 
the majority of mRNA species are classified as 



ranscript 



Putative change of expression following 
treatment according to dotblot 



Change according to RT-PCR 
quantitation 



Rat 



Guinea pig Rat 



Guinea pig 



lbumin 


N/A 


N/A 


No change 


No change 


ifunctional enzyme 


Up 


N/A 


Upregulated* (9 x ) 


N/O 


YP2C11 


Down 


N/A 


Downregulated* 


N/D 








(Abolished) 




YP4A1 


Up 


N/A 


Upregulated* (14 x ) 


N/D 


atalase 


N/A 


Up 


No change 


N/O 


D81 (TAPA-1) 


Up 


Down 


N/O 


Upregulated**(1.4 


ontrapsin-like protease inhibitor 


Down 


N/A 


Downregulated** 


x ) 
N/D 








(0.5 x) 




arathymosin-ot (zinc 2+ binding 


Down 


N/A 


Downregulated** 


N/D 


protein) 






(0.6 x ) 




ransferrin 


Up 


N/A 


Downregulated* 


No change 








(0.5 x) 




DP-Glucuronosyl transferase 


Down 


N/A 


Downregulated** 


N/O 








(0.2 x) 




ownUnknown-1 


Down 


N/A 


No change {P = 0.06) 


N/D 


n-ot2-glycoprotein 


Up 


N/A 


No change 


N/O 



u N/A. not applicable; N/O, not optimised; N/D, not done. 

*/><0.0005; 

** P<0.05. 



* 

v 

•'rare' in abundance (Bertioli et al., 1995). How- 

• ever, in their original paper describing the SSH 
technique, Gurskaya et al. (1996) demonstrated 
that SSH can enrich rare molecules between 1000- 

• and 5000-fold in a single round of hybridisation. 
Unfortunately, due to high background smearing 
in our initial experiments (which hindered identifi- 
cation of single bands), we were compelled to 
reduce the primary hybridisation time to only 4 h 
— a step that theoretically is likely to reduce the 
number of rare sequences (CLONTECHniques, 
1996). Furthermore, it has been claimed by the 
manufacturers that, whilst this technique can 
identify changes as small as 1.5-fold between the 
driver and tester populations, it is best suited to 
the isolation of genes that show a greater than 
5-fold increase (CLONTECHniques, 1996). In ad- 
dition, where tester and driver contain genes with 
large and small differences in abundance, the SSH 
method will be biased towards identifying those 
genes with the large differences (CLONTECH- 
niques, 1996). Thus, it is most probable that we 
have not isolated all of the more rarely expressed 
transcripts and those demonstrating small changes 

in expression. 

One problem that remains is identifying the 
function of genes isolated in SSH experiments as 
described herein, some of which may be crucial to 
the process of carcinogenesis, and are, to date, 
unidentified. However, we have provided evidence 
herein that SSH can be used to begin the process 
of characterising the extent and importance of 
altered gene expression in response to a chemical 
stimulus. The developments of this approach 
should include characterisation of temporal and 
dose responses, and functional analysis studies 
including knockout mice. In combination, such 
studies should make a significant contribution to 
our understanding of the molecular mechanisms 
of action and physiological relevance of gene reg- 
ulation in non-genotoxic hepatocarcinogenesis. It 
should then be possible to ascertain whether dif- 
ferentially expressed genes are causally or casually 
related to the chemical-induced toxicity, and 
therefore a substantial mechanistic advance. 

It is clear that there are also broader applica- 
tions for this experimental approach that go be- 
yond understanding the molecular mechanisms of 
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peroxisome-proliferator induced non-genotoxic 
hepatocarcinogenesis in rodents. The potential 
medical and therapeutic benefits of elucidating the 
molecular changes that occur in any given cell in 
progressing from the normal 40 the carcinogenic 
(or other diseased, abnormal or developmental) 
state are very substantial. Notwithstanding the 
lack of complete functional identification of al- 
tered gene expression, such gene profiling studies 
described herein essentially provides a 'fingerprint' 
of each stage of carcinogenesis, and should help in 
the elucidation of specific and sensitive biomark- 
ers for different types of cancer. Amongst other 
benefits, such fingerprints and biomarkers could 
help uncover differences in histologically identical 
cancers, and provide diagnostic tests for the earli- 
est stages of neoplasia. In addition, the genes 
identified by this approach may be incorporated 
into gene-chip DNA-arrays, thus providing a 
standard genetic fingerprint for a particular toxin 
treatment in a particular species. Interrogation of 
these gene arrays for an unknown compound that 
has a similar pattern to the known reference 
chemical would then provide evidence that the 
unknown may have a toxicity profile similar to 
the 'standard' fingerprint, thereby serving as a 
mechanistically relevant platform for further de- 
tailed investigations. 
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